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Abstract 

We study extremal surfaces in de Sitter space in the Poincare slicing in the upper 
patch, anchored on spatial subregions at the future boundary T + , restricted to constant 
boundary Euclidean time slices (focussing on strip subregions). We find real extremal 
surfaces of minimal area as the boundaries of past lightcone wedges of the subregions 
in question: these are null surfaces with vanishing area. We also find complex extremal 
surfaces as complex extrema of the area functional, and the area is not always real¬ 
valued. In dSi the area is real. The area has structural resemblance with entanglement 
entropy in a dual CFT. There are parallels with analytic continuation from the Ryu- 
Takayanagi expressions for holographic entanglement entropy in AdS. We also discuss 
extremal surfaces in the dS black brane and the de Sitter “bluewall” studied previously. 
The dS '4 black brane complex surfaces exhibit a real finite cutoff-independent extensive 
piece. In the bluewall geometry, there are real surfaces that go from one asymptotic 
universe to the other through the Cauchy horizons. 
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1 Introduction 

de Sitter space is fascinating for many reasons, in particular for holographic explorations to¬ 
wards addressing questions of cosmology and time. In this regard, some versions of dS/CFT 
duality [1] [2j 13] associate to de Sitter space a dual Euclidean CFT on the future timclike infin¬ 
ity X + boundary (in the Poincare slicing). A concrete realization in the context of higher spin 
theories appears in (4j. Further work on dS/CFT appears in e.g. [5] E* 3 E] [9. TOj [HJ 12j fTH] . 

In AdS/CFT , there has been considerable interest in understanding information theoretic 
notions in terms of geometric quantities via holography, in particular stemming from the Ryu- 
Takayanagi prescription [HUB] (see mm for reviews) for calculating holographic entangle¬ 
ment entropy of a subsystem in the strongly coupled boundary field theory. This is the area 
of a bulk minimal surface (in Planck units) anchored at the subsystem interface and dipping 
inwards upto a certain maximal depth typically called the turning point. A different way to 
think about this appears in [18]. More generally, these are extremal surfaces |19j . In this light, 
one might speculate that the bulk subregion enclosed by the entangling surface and the bound¬ 
ary subsystem in some sense encodes bulk physics corresponding to that part of the boundary 
theory contained in the subsystem, although a detailed understanding of the hologram (and 
bulk locality) would seem more intricate. 

It is interesting to consider these questions in the context of de Sitter space and dS/CFT. 
Assuming there is translation invariance with respect to a boundary Euclidean time direction, 
imagine constructing a subregion on a Euclidean time slice of the future boundary I + . Tracing 
out the complement of this subregion would lead to some loss of information and thereby give 
some associated entropy, which one might attribute to the subregion being entangled with 
the complement. In the bulk, intuition from the Ryu-Takayanagi prescription in AdS/CFT 
suggests that we study extremal surfaces in de Sitter space (in the Poincare slicing) on a constant 
boundary Euclidean time slice, defined as anchored on the subregion on the future (spacelike) 
boundary and dipping inwards (he. in the bulk time direction, towards the past). We find 
(sec. 2) that the bulk extremization problem exhibits some crucial sign differences from the 
AdS case. Focussing first on real surfaces, there are correspondingly some technical differences 
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such as the absence of a natural turning point (where the surface stops dipping inward). For 
sufficiently symmetric subregions such as strips (with an axis of symmetry), extremal surfaces 
can be defined as the union of two half-extremal-surfaces joined continuously but with a sharp 
cusp. Upon requiring that we choose minimal area, the extremal surfaces become null surfaces 
with zero area. In fact these are simply the boundaries of the past lightcone wedges of the 
subregion in question (restricted to the boundary Euclidean time slice), and are thus analogous 
to the causal wedges associated with causal holographic information [20j (note that these bulk 
causal wedges and the corresponding causal holographic information in general do not coincide 
with the bulk entangling subregion, and entanglement entropy). This answer - restrictions 
of past lightcone wedges - is well-defined for arbitrary boundary subregions, even without 
sufficiently high symmetry, and gives vanishing area. These surfaces with vanishing area do 
not appear to have any connection to entanglement in dS/CFT. 

It is therefore interesting to look for other extrema, in particular complex saddle points 
of the extremization problem, motivated by considerations in dS/CFT. For instance, in the 
formulation of [3] of the dS/CFT dictionary, the CFT partition function is Zcft — V I ; where 
T is the bulk late time Hartlc-Hawking wavefunction of the universe subject to appropriate 
(Bunch-Davies) boundary conditions at early times. In a semiclassical approximation T ~ e lScl , 
the dual CFT energy-momentum tensor ( TT) correlators exhibit central charge coefficients of 

ftd-l 

the form Cd ~ i 1 ~ d G^ ( w ^ c ^ essentially an analytic continuation from Euclidean AdSj+i)- 
With this in mind, focussing again on strip subregions in the present context, we indeed find 
these complex extremal surfaces: they exhibit “turning points” in the interior. They should 
be thought of as living in some auxiliary space, and are distinct from the bulk past lightcone 
wedges (which define real subregions in bulk dS^). The area of these surfaces is in general not 
real-valued. In dS±, we find that x(r) parametrizing the strip width being real-valued suggests 
that the bulk time r parametrizes a complex path r = iT. The area (in Planck units) of these 
complex surfaces in dS± is real-valued and negative, while in dSd+i with d even, the nature of 
these surfaces is different and the area is pure imaginary. The area has structural resemblance 
with entanglement entropy in a dual (non-unitary) CFT^, with central charge Cd ~ i 1 ~ d -Q^ , 
with a leading area law divergence, and subleading terms. There are parallels with analytic 
continuation from the Ryu-Takayanagi holographic entanglement expressions from AdS. It 
is a useful consistency check that these central charges Cd ~ i 1 here resemble those in 

the (TT) correlators in [3], mentioned above. From the point of view of the dual Euclidean 
CFT living on the future boundary I + , one might formally associate a density matrix w.r.t. 
boundary Euclidean time evolution and a reduced density matrix to the subregion obtained 
by tracing out the complement. It would be interesting to explore this further, perhaps in 
dS/CFT as entanglement entropy in the dual Euclidean CFT. 

We then discuss (sec. 3) an asymptotically de Sitter space [13] - the dS black brane - 
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where subleading normalizable metric components are turned on: in dS^/CFT 3 , they have the 
interpretation of saddle points representing the Euclidean CFT with uniform energy-momentum 
density expectation value. The corresponding extremal surfaces in the d,S± black brane exhibit 
a finite cutoff-independent real-valued extensive piece (again negative) with some resemblance 
to a thermal entropy. Finally we discuss (real) extremal surfaces in the closely related dS 
“bluewair geometry, which are not obtained by analytic continuation: there are real extremal 
surfaces which cross from one asymptotic universe to the other through the Cauchy horizons. 

2 Extremal surfaces in de Sitter space 

de Sitter space dS d+ 1 in the Poincare slicing or planar coordinate foliation is given by the metric 

ds 2 = ^% L (— dr 2 + dw 2 + dx 2 ) , (1) 

where half of the spacetime, e.g. the upper patch, has Z + at r = 0 and a coordinate horizon at 
r = — oo. This may be obtained by analytic continuation of a Poincare slicing of AdS, 

r ->• -it , R AdS ->• -iRds , t -> -hr , (2) 

where w is akin to boundary Euclidean time, continued from time in AdS. 

The dual Euclidean CFT is taken as living on the future r = 0 boundary X + . We assume 
translation invariance with respect to a boundary Euclidean time direction, say w, and consider 
a subregion on a w — const slice of Z + . One might imagine that tracing out the complement of 
this subregion then gives entropy in some sense stemming from the information lost. In the bulk, 
we study de Sitter extremal surfaces on the w = const slice, analogous to the Ryu-Takayanagi 
prescription in AdS /CFT. Operationally these extremal surfaces begin at the interface of the 
subsystem (or subregion) and dip inwards (towards the past, in the bulk time direction). For 
simplicity, consider a strip on the w = const surface (i.e. a constant boundary Euclidean time 
surface): this bulk d-dirn subspace has metric 

ds 2 = ( - dr 2 + dxfj • (3) 

Xi^w 

This is not a spacelike subspace in the bulk and it might seem that the extremal surfaces are 
timelike in general: however we will find that this is not the case. 

2.1 Real extremal surfaces 

Let us consider a strip subregion with width direction say x, the remaining Xi being labelled 
y l . A bulk surface on the w = const slice bounding this subregion and dipping inward (towards 
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the past) is bulk codim-2: its area functional in Planck units is 


S ds = 


4 G. 


d +1 


, d—2 

n 

i=l 


Rdsdyi Rds 


\fd 
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dx 2 
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Rjs^ti 

4Grf + i 



/'dx\ 2 

\dr) 


( 4 ) 


We consider extremizing the action to find extremal surfaces with minimal area, along the lines 
of the Ryu-Takayanagi prescription for entanglement entropy in AdS. The S d s extremization 
gives a conserved quantity (x = jp) 


x 

\/l — x 2 


Br d ~ l 


x 


1 

1 + £?2 r 2d-2 


i.e. 


x 


j^2 T 2d-2 

1 + B 2 r 2d ~ 2 


(5) 


We see that x 2 —> 0 near the boundary r —>■ 0. Assuming the conserved constant satisfies 
B 2 > 0 makes all the expressions real-valued and means x 2 > 0, with x 2 —» 1 in the deep 
interior for large |r|. For B 2 > 0, these are timelike surface^]. This gives the solution (upto 
boundary conditions) and corresponding area integral 


x{t) = ± 


Br d ~ x dT 
yj\ + B 2 T 2d ~^ 


±X(r) , 


Rts 1 t / r _ 

4G d+ i d 2 /. T- 1 ' 1 v 2 ! + B 2 r 2d ~ 2 ' 


( 6 ) 


The main difference between this case and the minimal surface in AdS stems from B 2 > 0 
implying that there is no smooth “turning point” where x 2 = t co. In fact B 2 > 0 

means x 2 is bounded, with 0 < x 2 < 1. For any finite B 2 > 0, the extremal surface in this case 
begins to dip inwards from one boundary of the strip subregion and (rather than turning around 
as in AdS ) continues indefinitely, eventually approaching x —> ±1. With a view to associating 
a bulk subregion with the boundary subregion in question, let us artifically cut off the inward 
dipping surface at some interior location r = r 0 , the bulk subregion then defined by the interior 
of the boundary strip subregion and the joined surface. So consider the half-extremal-surfaces, 


x l (t) = X(t) - X(t 0 ) = - x r (t ) , 

xl(tq) = 0 = x r (t 0 ), x l ( 0) = - l - = - Xfl ( 0 ) =► l - = X(t 0 ) . ( 7 ) 

This gives an extremal surface made of two half-extremal-surfaces joined continuously but with 
a sharp cusp at To (see Figure ED- This defines the corresponding wedge-like bulk subregion, 
enclosed by this extremal surface and the boundary subregion. These conditions do not de¬ 
termine the parameters B, To uniquely, given the subregion width l. Varying B gives different 
extremal surfaces. By comparison, in the AdS case, the turning point r* = is fixed by the 
global nature of the entangling surface as the location where x 2 — > oo, the surface turning 
around. 

1 One might instead want to consider spacelike surfaces with x 2 > 1 and therefore take, instead of dU), the 

area functional as Sds = I \j ( 37) 2 ~ 1- We will discuss this in the next subsection. 
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To follow the Ryu-Takayanagi prescription, we would want to identify those extremal sur¬ 
faces that have minimal area@. From ()6j), we see that as B increases, the area S d s decreases. 
Furthermore, (JSJ) shows that as B increases, x 2 increases and eventually approaches x 2 —> 1 as 
B —> oo. In this limit, x{r) —>■ ±r and Sds appears to vanish. In fact this is a sensible result: 
in hindsight, it should have been obvious from ([4)) that minimal area arises when the extremal 
surface becomes null. This null extremal surface is in fact simply the boundary of the past bulk 
lightcone of the subregion, restricted to the boundary Euclidean time slice. 



Figure 1: Extremal surfaces in de Sitter made of two half-extremal surfaces joined continuously but with a 
sharp cusp at r 0 . As B increases (till eventually B ^hr)> the surface approaches x 2 —> 1 (figure on the right). 


An alternative argument corroborating the above conclusion is the following. Physically, 
the shortest length (or time) scale here is Tuv = e so that in ([5]) when Be^ 1 1 we can 


approximate x 2 ~ 1 and so x(r) ~ ±r giving | ~ tq. Thus one might estimate (EJ) as 


S, 


dS 


pa!—1 

dS Vd -2 
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dr 


4 G, 


d +1 
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d -1 
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4 G, 
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T d ~i Vl + B 2 r 2d ~ 2 
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( 8 ) 


where the second integral can be seen to vanish as r —> 0. Now the first integral scales 
as while the second integral can be expressed in terms of the hypergeometric function 


£>d— 1 

2 -Fi as -42-■ 


V d -2 


r[l - 2 FiG 


d—2 


-B 2 T 2d ~ 2 


Gd+1 {d -2)r ^ L " ““ _ 2 (d 3 i)’ '0 )] ( the extremal surface in © is itself 

expressed as x(r) = ±r 2 ^i(|, j. im- §5^, ~ or ±Br d 2 Ti(|, |^|, -B 2 T 2d ~ 2 ), using 

the integral representations of 2 Fi). As B 2 increases, this second integral is seen to scale as 

R'' 


^V d _ 2 B {d - 2)/ ( d - 1) . Thus when B 
G d + i a 2 


3^1 this cancels the earlier contribution and we again 
see the leading S d s scaling to be vanishing. 

In dS 3 (i.e. d = 2), we obtain x(r) — ± f ^^ d l T , = ±-^\/l + B 2 r 2 and the boundary 
conditions give ^(y/l + B 2 Tq — 1) = |, the area integral becoming S d s = X To y 1 n/1+ ] b2t2 • 
Analysing these vindicates the conclusions above. 

It is worth noting that our construction of joining two half-extremal-surfaces appears invalid 
unless the subsystem has sufficiently high symmetry (in particular an axis of symmetry). Rclat- 
edly, one might look askance at the entire extremization procedure here, in particular whether 


2 Note that surfaces with maximum area correspond to minimizing B: this gives B = 0, which are disconnected 

R ‘cis 1 y r° dT 

4 G d+1 v d—2 J £ 


surfaces x(t) = const , with area Sds 


A with a leading divergent piece. 


5 
































one allows non-smooth surfaces (with cusps) in the extremization: one might then wonder if 
more general surfaces need to be considered, e.g. a zigzag null surface formed by joining multi¬ 
ple partial surfaces with multiple cusps. This would be useful to systematise more rigorously. 
However the final answer, the past lightcone wedge, is well-defined for an arbitrary subregion, 
comprising two piecewise smooth extremal surfaces joined with just a single cusp (rather than 
multiple cusps). The past lightcone wedge boundary (restricted to the boundary Euclidean 
time slice) is however a complicated surface: it would be interesting to understand the shape 
dependence here. The resulting area is of course always zero for all these null surfaces, and 
does not reflect entanglement structure. 

From the point of view of bulk de Sitter alone, one could consider volume subregions in 
the full d-dirn boundary Z + (at r = 0), i.e. not on the constant boundary Euclidean time 
slice. These would give codim-1 surfaces. For a strip subregion with width direction say x, the 
remaining X{ being labelled i/i, analysing the area integral of the bulk surface for extremization 
gives 


S. 


dS 


dr 




/dxy 
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2^2 d 


x 2 = 


B 2 T 


1 + B 2 T 2d ’ 


Sds = RdsVd-i 


VT 

r T o 


x 2 


dr 


= Br d 
1 


T d y/1 + B 2 r 2d ' 


(9) 


This has volume scaling. Again S^s decreases with increasing B, with x 2 —» 1: the resulting 
extremal surfaces are null surfaces defining the past lightcone wedges of the volume subregion, 
with vanishing area. 

These real null surfaces with vanishing area do not appear to have any bearing on entan¬ 
glement in dS/CFT. In what follows, we will explore other complex saddle points. 


2.2 Complex extremal surfaces 

For what follows, it is useful to recall the dS/CFT correspondence, for de Sitter space in 
the Poincare slicing (JT]) , obtained by analytic continuation ([2]) of Poincare AdS. A version of 
dS/CFT [H El Sj states that quantum gravity in de Sitter space is dual to a Euclidean CFT 
living on the future boundary Z + . More specifically, the CFT partition function with specified 
sources 4>io(x) coupled to operators Oi is identified [3] with the bulk Hartle-Hawking “wave- 
function of the universe” as a functional of the boundary values of the fields dual to Oi given 
by 4>io(x ). In a semiclassical approximation, this becomes Zcft = v h[0jo(^)] ~ e lScl ^ i0 ^ where 
we need to impose regularity conditions on the past cosmological horizon r —> — oo: e.g. scalar 
modes satisfy 4>k( r ) ~ e lkT , which are Hartle-Hawking (or Bunch-Davies) initial conditions. 
Operationally, certain dS/CFT observables can be obtained by analytic continuation (J2]) from 
AdS (see e.g. [3j, as well as 0)- The Bunch-Davies initial condition itself can be thought 
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of as analytic continuation of regularity in the AdS interior. Dual CFT correlation functions 
can be obtained from the dictionary Zqft — '&■ For instance (from [3]) in a semiclassical 
approximation T ~ e lScl , a massless scalar in dS 4 has a mode solution 0 = {i~ikr)e ik ^ w ^h 
Bunch-Davies initial conditions </> ~ e tkr at early times (|r| —>■ oo) and r c a late-time cutoff. 
The classical action evaluated on this solution is iSds 4 ~ J d 3 k(iy- — k 3 + .. .) 0 °-dji , the 
divergent terms being oscillatory (pure imaginary). Appropriate graviton modes can be approx¬ 
imated as massless scalars so that dual CFT energy-momentum tensor ( TT) correlators can be 
read off as ( T(k)T(k'))ds 4 = 0^p"l</>°=o ~ (~%f)k 3 8 (k + k'). The central charge coefficient 

here is C 4 ~ essentially an analytic continuation from Euclidean AdS 4 . I 11 dS$ [3], we have 
iSds s f d 4 k(f ) 0 ■ - + k A log(-r c fc) + ...) where we have only exhibited the nonlocal term 

which contributes to the 2 -point function: this gives ( T{k)T(k , ))ds 5 ~ (iyrr)fc 4 log k 8 {k + k'), 

3 —1 

with central charge coefficient C 5 ~ ^~gT- More generally we have Cd ~ i 1 d G dS (which is 
essentially an analytic continuation from EAdSd+i)- 

From this dS/CFT point of view, one might expect any leading divergence in the CFT 
entanglement entropy (assuming it exists) to be of the form Cd^Ar which is in general complex¬ 
valued: thus it is natural to ask if there are additional (perhaps complejj§) extrema of the area 
functional that should be considered in de Sitter space, with possible dS/CFT interpretations. 
With a view to considering spacelike surfaces with x 1 > 1, let us take, instead of ([4j), the dSd+i 
area functional on a w = const slice as 


Sds = 


4Grf + i 


dr I/dx \ 2 

T d -i VUJ 


-1 , 


x 

sjx 1 — 1 


At 


d -1 


( 10 ) 


where we are considering strip subsystems with width along x. The second expression above 
is the conserved quantity obtained in the extremization. This is essentially the same as ([5]), 
but with B 2 = —A 2 < 0, and A being real. One might ask if this can be interpreted as a 
real surface x 2 = ^- 53 = 23 r • However we need to require that the surface reaches the boundary 
r —» 0 from where it drops down (inward): near r —> 0, we find x 2 rs-/ —A 2 T 2d 2 , so that 


_^2 r 2d-2 

2 _ j^2j2d—2 ’ 


(ll) 


this being a complex surface in some sense. 

Let us focus now on dS^/CFT^, for concreteness, to understand this better. The extremal 
surface near r —>■ 0 in this case is x{r) ~ ±iAr 3 + x(0). We want x(r) to be real-valued since it 
parametrizes a space direction in the dual CFT 3 : this requires that r takes imaginary values. 
In more detail, near r —» 0, we have x — > ±| and the two ends of the surface are parametrized 

3 Recall that complex geodesics appeared in [2T] in the context of the black hole interior. Complex extremal 

surfaces have recently appeared in a different context [22], as well as [23] ■ 


7 














rsj 


as xl(t) ~ | + iAr 3 and xr(t) ~ | — iAr 3 . For to be real-valued with A real, we 

must have pure imaginary r = iT with T real, giving xl ~ | + AT 3 ~ — as T increases, 

Xl increases from — | and Xij decreases from The global structure of the surface shows a 
“turning point” at t*A 2 = 1, where x 2 —> oo, very similar to the situation in AdSl From the 
point of view of the discussion in the previous subsection, the two half-extremal-surfaces xl,xr 
in this case join smoothly at the turning point r* as in AdS, with Xl(t*) = 0 = x#(t*) and 
Xl, xr matching. This gives the width 


Ax 

~Y 
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At 2 dr 
v/1 — A 2 t 4 


t=iT 


" T * i.i 3 .AT 2 dT 
) Vl - A 2 T 4 


= =>■ r* ~ il . 


( 12 ) 


The reality of Ax = l with A real again suggests that we parametrize the r-integral over the 
path t = iT in a complex r-plant^: we have then rescaled T using A to make the integration 
variable dimensionless (and # = / 0 ^= = -v/TTpry )• The turning point here is t* = -j=. The 


yi-y 


VA- 


integral can be parametrized in terms of hypergeometric functions 2 -Fi. The extremal surface 
x(t) with r imaginary does not correspond to any real bulk subregion in dS± enclosed by 
the surface, but really lives in some auxiliary space. In a sense, the structure here is very 
much like analytic continuation of the AdS 4 expressions a la Ryu-Takayanagi: we will discuss 
this more below. From that point of view, since the analytic continuation (J2]) faithfully maps 
AdS 4 -H- dS 4 , this is a faithful map from the subsystem to the auxiliary bulk subregion. The 
area now becomes 


S df! = 2 ^kv, 
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.R 2 * 


4G d 
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T 2 y'lV 




dS 
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dr 
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In principle, we could assign ±i in the second expression, as a choice of the branch of the 
square root: the choice of the minus sign leads to an appropriate coefficient as we see below. 
The integral itself is just as in AdS 4 , giving 


& = 


p / 3 \ 2 

where c 3 = 27r(j^fy) is a constant as in AdS, stemming the finite cutoff-independent part of 
the integral. Note that here we have used the relation tuv — R f° r the ultraviolet cutoff in 


1 1 

-c 3 —] = 

Tuv T * 


R dS y 

2 G 4 li U 3 l 


CV(I-c 3 i 


(14) 


4 Strictly speaking, this may be too restrictive. We have required a;(r) for all r be real-valued: this means 

that each point on the surface directly maps to a corresponding real-valued spatial location within the strip in 

the dual CFT 3 . One might instead think that one need only require the boundary value x(0) be real, which 

would not restrict the r-path. This would suggest more general complex extremal surfaces defined over complex 

r-space, with the width Aa: required to be real-valued. See e.g. |23j for some discussions along these lines: I 

thank S. Fischetti for a discussion on this point. 






















the dual Euclidean field theory, suggested by previous investigation^ in dS/CFT. Also we 
have rewritten the last expression in (1141) in terms of the dS 4 /C FT :i central charge C 3 ~ ~~cf 
appearing in the (TT) correlators in [3], reviewed above. Sds hi (TT4|) is real-valued and bears 
structural resemblance to entanglement entropy in a dual CFT 3 with central charge C 3 ~ 
~~cf < The first term resembles an area law divergence [321 33], proportional to the area of 
the interface between the subregion and the environment, in units of the ultraviolet cutoff. It 
is also proportional to the central charge which represents the number of degrees of freedom in 
the dual CFT: in this case, C < 0 reflecting the fact that the CFT is non-unitary. The second 
term is a finite cutoff independent piece. Whether the expression (UTD physically is holographic 
entanglement entropy in dS^/CFTs is less clear from this bulk extremal surface analysis. 

In some sense, —S d s appears to resemble entanglement entropy in AdS/CFT , sharing var¬ 
ious features including subadditivity. For instance, the quantity I d s[A, B] = S d s [A] + — 

S ( is [A U B] for two disjoint subsystems A, B , exhibits various properties of holographic mutual 
information in AdS including an analog of the disentangling transition in the classical gravity 
approximation [34], but with some crucial differences. For strip subregions that are sufficiently 
nearby but disjoint, I d s[A,B\ is nonzero: e.g. using (U4|) for a single strip, we obtain for two 
parallel strips of equal width l and separation x, 


I<is[A, B] — S dS [A] + Sdsf-B] — S dS [A U B] 


Rls T W 2 1 1 

--p^—C 3 Vi ( — y -b 

G 4 V l x 21 + x 


(15) 


U B] arises from the area of the connected surface between A, B as S(2l + x) + S(x). 
This is similar to the structure of holographic mutual information for strips in AdS 4 , e.g. 
the UV divergent pieces cancel, with a cutoff-independent divergence C -y 1 as the subregions 
collide. The striking difference is that I d s[A,B] < 0, rather than positive definite, following 
from the fact that C = < 0. Thus I d s[A, B] is large and negative when the subregions are 

nearby, then increases as the separation | increases, and eventually approaches zero as | —> 
v ^~ 1 ~ 0.62. Beyond this critical value, I d s[A,B] vanishes identically and the two subregions 
are disconnected. This disentangling transition in the classical gravity approximation arises 
from the transition between the connected and disconnected surface for A U B. What we are 
seeing is that Sds[A] -bSdsf-B] < S d s[AU B], i.e. —S d s satisfies strong subadditivity for disjoint 
parallel strips A, B. By comparison, using the real lightcone wedge surfaces in the previous 

5 See e.g. glia sun], which discuss this (in some cases implicitly). Heuristically, we expect that evolution in 
the bulk direction is encoded by renormalization group flow in the dual field theory: see e.g. [24 f (25l 1261 l27l f28f i29 ] 
and more recently e.g. [Ml El] for discussions on this in the AdS context. In the present dS case, the bulk 
description is time evolution i P- = 'H'fr, with H. being an evolution operator. Through dS/CFT, this becomes 
i = T-LZcft- Heuristically this maps to a renormalization group equation schematically of the form 
5 z °f t = OZqft, if Tuv = *£) with O an appropriate operator generating RG flow. If we view e as a floating 
RG parameter, this again suggests the path r = iT in complex r-space for a dS/CFT interpretation. 
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subsection, we see that disjoint boundary subregions are always causally disconnected and thus 
uncorrelated for any nonzero separation. Correlation functions are nonzero: the disentangling 
transition above is in the classical gravity approximation, and we expect subleading terms in a 
large-N-like expansion of I d s[A,B] (see e.g. [35j in the AdS context). 

We now discuss dS d +i for even d (in particular dS 3 , dS 3 ) where the nature of these extremal 
surfaces is different. We would like to retain the relation tjjv — as following quite generally 
in dS/CFT from time evolution mapping to renormalization group flow. This suggests we 
parametrize the bulk time parameter r along a complex path r = iT as in dS±. However now 
with A 2 > 0 the surface (HD near r —> 0 gives x ~ ±iAr d 1 , i.e. x(t) ~ ±iAr d . Thus x 
cannot be made real-valued for any even d in this manner. A way out is to take the parameter 
A 2 —y —A 2 : the surface equation now is the same as (J5]) but with the bulk time parametrized 
as r = iT. The expressions (ITU . (TlTh then give 


j^ 2 r 2 d -2 

1 + A 2 r 2rf-2 ’ 


S ds 


. Rgs 1 v [ r * dr 2 

*4G d+ i d ~ 2 J TUV r d ~ l Vl + A 2 r 2d ~ 2 

^ s-^v d 2 r— _-_ 

2 G d+ i J e T d ~ l + (_iy-i A 2 T 2d - 2 


. (16) 


For even d, the (—l) d 1 gives rise to a “turning point” at T 2d 2 A 2 = 1: the width now scales 
as l ~ T* ~ — fr*. The integral is as in AdS, giving 


S ds 


r^i 


•1 -d 


R' 


id—1 

dS 


2 G, 


-V d -2 


d+1 


1 

e d-2 



(17) 


The leading divergence S' 


div 

dS 


~ l 


r>d — 1 

•1—dfCs 
1 G d+ 


jE^d — 1 t r 

i l ~ d 7rSF — 3=7 resembles an area law: the central charges C d 

^^d +1 e 


- here resemble those in the ( TT) correlators in [3j reviewed above. This leading be¬ 
haviour appears independent of the shape of the subregion, expanding (fTOl) and assuming that 
x is small near the boundary tuv■ Unlike dS^, note that S d s in dS d +1 with even d is not 
real-valued, in particular for dS 3 ,dS 5 . For instance, in dS 3 , we obtain from (fIBD 


r = iT , x ( t ) ~ ±^-VTTkUr 2 , S ds ~ log — = log - . (18) 

A G 3 Tuv G 3 e 

Note that x{t) appears real, although the parametrization is r = iT. 

It is interesting to recall the Ryu-Takayanagi expression for entanglement entropy for an 
(infinitely long) strip-shaped subsystem with width along the x-direction, given as the area of 
the corresponding minimal surface in the bulk AdS d +1 geometry (with radius R), 


S A ds[R,x(r),r} 


4G*d + i 


Ud-2 



/dx \ 2 

\dr) 


j \2 


(x 


j^ 2 r 2 d -2 
^ _ J^2y2d—2 1 


(19) 


10 























where the conserved quantity A in the extremization is related to the turning point as r f -1 = A. 
Noting that dS d +i in Poincare slicing (J3]) is just the analytic continuation of the corresponding 
t = constant spatial slice in AdS d +i, obtained by (j2J), be. r — y —ir, t — > —iw, R —> —iR d s, let 
us carry out this analytic continuation on the Ryu-Takayanagi expression. Indeed we see that 
Sds in (J4]) appears very much like the analytic continuation of SAds[x(r),r], with the various 
factors of i conspiring to leave a single i behind, be. 


S AdS [R,x(r),r\ ->■ -i 


R 


d— 1 
dS 


-Vd- 


d—2 


= S d s{Rds,x(r),T} 


4G d+ i J t“ 

On the analytic continuation of the extremization itself, we obtain 

( _ ■\2d—2j^2,j.2d—2 _ / ■y'jd-l Jp r 2d-2 


( 20 ) 


— X 2 = 


l _ (j\2d-2j^2 T 2d-2 ’ 


i.e. 


X 2 = 


1 — ( — 1 Y~lj\2 T 2d-2 ’ 


SdS — 


R 


•d -1 
dS 


4G, 


-Vb- 


dr 


d—2 


d +1 


T' 


d -1 


V'l - (-1 y-iA 2 T 2d - 2 ' 


( 21 ) 


This expression corroborates the minus sign in (TT 6 l) and (1T3]) . (TT4l) . The analytic continuation 
essentially recovers our earlier calculations in dS^ and dSd+i for even d. For instance, in dS$ 
(i.e. d = 4), we obtain 


A 2 r 6 
1 + A 2 t 6 


S dS = 


i^V 2 

4G, 


dr 1 

r 3 VI + 


( 22 ) 


With real A, this is as such a real extremal surface as in the previous subsection: taking A large 
minimises the area and we obtain the null surfaces earlier with vanishing area representing the 
past lightcone wedge of the subregion. However parametrizing as r = zT, there is a turning 
point at r* = ^ 73 , and a corresponding complex surface and corresponding area given by (TT 6 |) . 
The area in (l 22 ]l then becomes S d s ~ The extra i can be thought of as arising 

from the odd powers of Rds under the analytic continuation from AdS 5 . Thus interestingly for 
even d (in particular, dSs and dSs), the expression S d s obtained by analytic continuation of the 
Ryu-Takayanagi entanglement prescription leads to complex surfaces with corresponding area 
Sds pure imaginary. 

To summarize, we have studied bulk de Sitter codim-2 extremal surfaces. Real extremal 
surfaces are the boundaries of the past lightcone wedges of the boundary subregions, with 
vanishing area. Codim-2 complex extremal surfaces have area exhibiting structural resem- 

R d ~ i n 

blance with entanglement entropy in a dual CFT^ with C d ~ i 1 ~ d 7 T S — central chargqj matching 
those appearing in the (TT) correlators obtained from the wavefunction of the universe [3]. 

6 Note that codim-1 surfaces with area functional S ~ Rds l \ d ^ 1 f fjVx 2 — 1 upon extremization exhibit area 
whose leading divergence does not reflect the CFT central charge. 
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In (IS 4 / CFT :i , the area is real-valued and negative: in this sense, these complex surfaces have 
lower area, suggesting that they are the preferred minimal surfaces. Our calculations here have 
been done for a strip subregion but it would appear that generalizations to other subregion 
shapes will exhibit similar features. For instance, the spherical subregion extremal surface pre¬ 
sumably exhibits a logarithmic term with interesting coefficient (this universal coefficient in the 
logarithmic term has been studied recently in [36], exhibiting agreement with the correspond¬ 
ing coefficient in the logarithmically divergent term in the wavefunction of the universe in the 
classical approximation). 

It is worth noting that this analysis of bulk extremal surfaces is different from studies of 
entanglement entropy of bulk fields in de Sitter space e.g. [37, 381 139, 3D]- 


3 Extremal surfaces in the dS black brane 


We now study extremal surfaces in the asymptotically dS spacetime studied in [13], i.e. 


ds 2 



dr 2 

1 + ar$T d 


d -1 

+ (1 + ar d r d )dw 2 + ^2 d' x t) j 

i =1 


(23) 


with a a complex phase and 7"o is some real parameter of dimension length inverse. An analog 
of regularity in the interior for an asymptotically AdS solution is obtained here by a Wick 
rotation r — il and demanding that the resulting spacetime (thought of as a saddle point 
in a path integral) in the interior approaches flat Euclidean space in the (/, ra)-plane with no 
conical singularity. This makes the w-coordinate angular with fixed periodicity (and l is a 
radial coordinate), giving a = l > r 0 , w ~ w + Thus the spacetime (1231) 

is a complex metric which satisfies Einstein’s equation with a positive cosmological constant 
Rmn = S~gMN , A = • This resulting metric satisfying regularity is equivalent to one 

obtained by analytically continuing the Euclidean AdS black brane 


ds 2 


R Ads f dr2 
r 2 V1 - r$r d 


d -1 


+ (1 ^ r d r d )d9 2 + J2 dx l), 

i=l 


(24) 


where 6 ~ 6 + ^i)r 0 > ^ ie asymptotically de Sitter spacetime (123|) using (]2]) and we identify 
r 0 = To, giving the phase jzjyi- Th e regularity criterion is simply the analog of regularity of 
the EAdS black brane. The condition l > t 0 is equivalent to the radial coordinate having the 
range r > ro- We see that “normalizable” metric pieces are turned on in (1231) . We then expect 
a nonzero expectation value for the energy-momentum tensor here, as in the AdS context 
pEl , HU 12E, l2Hj- In the present case [T3], we have T tJ = 5Z s %fj T = cc 9if ■> w l iere 

gdj* is the coefficient of the normalizable r d ~ 2 term in the Fefferman-Graham expansion of the 
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metric (123|) . This definition of T tJ is natural for a CFT with partition function Zcft , equated 
with T: thus, most notably, the i arising from T, the wavefunction of the universe, implies 
that the energy-momentum tensor is real only if gff is pure imaginary. In effect, this dS/CFT 
energy-momentum tensor can be thought of as the analytic continuation of the EAdS one. The 
spacetime (T23T) for dS 4 /CFT 3 gives real T^, with T ww = — with T ww + {d — l)T tl = 0. 

The w-coordinate is naturally interpreted as Euclidean time from the structure of the energy- 
momentum tensor: so let us now consider a strip subregion on aw — const surface in (123|) . The 
area functional (in Planck units) of a bulk surface bounding this strip and dipping inwards is 


S ds 


; R dS l y . [ I 1 ( dx V 

Gd+1 d 2 J r d_1 y 1 + ctTy T d \dr) 


(25) 


defined so that for r 0 = 0, this reduces to our de Sitter discussion in sec. 2.2. For the dS 4 brane 
(i.e. d — 3), we obtain for the extremization, 


^ = r iA ^ 2 dr s = _yj§s r dj_ _ 2 _ ^ 

Jo yj(l — r 3 ) (l — A 2 t 4 ) ^ 4 Jwv T \J (l — ir^r 3 ) (l — A 2 r 4 ) 

Now for small width /, this is essentially similar to the previous discussion on pure dS 4 and 
we have ^ ^ ~ r* = ^=, where A is real. In particular, the width Ax being real-valued 

suggests that r parametrizes a complex path r = iT with T real. As l increases however, 
the other denominator approaches a zero also, with r —> In this limit, we thus have 
r* —* ~ ^ and large l ~ — ir*, obtained from the double zero as 


A X = l_~ r* iAt ? dr r=l T f T * i.i^ATl dT 

2 2 J 0 — ir$T 3 )( 1 — A 2 t a ) Jo a/(1 — 7 oT 3 )(1 — A 2 T 4 ) 

Note that reality of the width Ax implies now that the range of T is restricted as T < A he. 
asymptotically r —> This is similar to the fact that in the AdS black brane, static minimal 
surfaces in the IR limit (large subsystem width) wrap the horizon but do not penetrate beyond. 

Now the area integral exhibits a cutoff-independent piece which can be estimated from the 
contribution in the deep interior where r —> r*: the contribution to the integral near the double 
zero thus scales as iAx giving 


S fin ~ 


-Mk 1 

G 4 T -2 


(<0 


d 2 
n dS 


g 4 


Tq VlI = 


— f r r 2 


CT‘V,I 


(28) 


j ^2 

which resembles an extensive thermal entropy in a 3-dim CFT with central charge C ~ 
at temperature T 0 = r 0 . Note that S An < 0. In fact S An is the analytic continuation ([2]) of 
the horizon entropy -^- TqViI of the EAdS 4 black brane (|24|) . which can be obtained as the 
horizon contribution from the partition function in the classical approximation. 
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We recall that the entanglement entropy area functional for the AdSj+i black brane from 

Vd - 2 - R d_1 r dr 


(29) 


the Ryu-Takayanagi prescription is S' = f fa d r x ) 2 + jfaj, giving 

, i.‘) A 2 r 2d ~ 2 „ V^ 2 R d ~ 1 r* dr 2 

{x ') 2 = 


s = 


(1 - r^r d )(l - A 2 r 2d ~ 2 y ~ 4G d+ i J e r d ~ l ^/(i - r ((r d )(l - A 2 r 2d ~ 2 )' 
Under the analytic continuation, we obtain 

A 2 T 2d ~ 2 


X 2 = 




(l — (— i) d TQ T d ) (l — ( —l) d 1 A 2 T 2d 2 ) ’ 
d-2R d dS l r ~ idT 2 


4Gd+l J Tuv T \J~ (—i) d T$T d ) (l — ( — l) d_1 H 2 T 2d_2 ) 

For generic dimension d, we see that S is not real, as in the earlier discussion with To = 0. 


(30) 


3.1 The de Sitter bluewall 

We now explore metrics of the form 


, but with the parameter a = — 1 


herc0, i.e. 




dr 2 


1 s-r-dsT-d 

1 /q / 


+ (1 - T d T d )dw 2 + dx \= 


R 


■dS 


dr 2 

7F: 


+ f{r)dw 2 + dx\ 


(31) 


The w-coordinate here has the range —oo < w < oo. This spacetime [13] has a Penrose diagram 
shown in Figure [2] which resembles that of the AdS black brane rotated by there are two 
asymptotically dS universes (for r < 4-), and timelike singularities cloaked by Cauchy horizons 
at r = 4., which “cross” at a bifurcation region. The Penrose diagram has many similarities 
with the interior of the Reissner-Nordstrom black hole (or wormhole). Late time infalling ob¬ 
servers near the Cauchy horizon see incoming lightrays from early times as highly blueshifted, 
essentially stemming from lightrays “crowding” near the Cauchy horizon, suggesting an instabil¬ 
ity. It is unclear if this spacetime has any interpretation in dS/CFT : nevertheless, formally, one 
finds the energy-momentum T tJ to be imaginary in dS^/CFT^, perhaps reflecting the blueshift 
instability here. Here we will simply look for bulk codim-2 extremal surfaces, lying either on a 
w = const slice or an x = const slice (from a bulk point of view alone, either could be taken as 
Euclidean time slices), restricting to real surfaces which may also be timelike. 

The area functional for a surface in (IHTj) bounding a subregion on a x = const slice of Z + , 
and wrapping the other x t ^ x, is 


S = 


pd—i 


4 G 


■<2+l 


dr 

-t-cL—1 




.d^d 


f M = 1 - tJt 


(32) 


7 The metric (l23l) with a = +(—i) d is similar to the dS black brane, except with T t j of the opposite sign, 
while a = +1 gives a real spacetime with a spacelike singularity at 


■ oo. 
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This does not correspond to any analytic continuation from the Ryu-Takayanagi formula for 
the AdS 4 black brane, so we analyse this directly focussing on the dS 4 bluewall. Along our 
earlier discussions in sec. 2, we find real extremal surfaces corresponding to 

^ _ 1 f dT 1 (' 33 ') 

(1 - TqT 3 ) 2 1 - TgT 3 + B 2 t 4 ’ 4 G 4 J T d_1 - r 3 r 3 + B 2 r 4 ’ 

where the constant R arises from a conserved quantity in the extremization. The first equation 
describing the surface can be rewritten as 

_ dr (dw\ 2 _ , 2 B 2 t 4 

1 — Tq r 3 ’ Vdr*/ * 1 — TqT 3 + R 2 r 4 

where we are using r* here for the “tortoise” coordinate in this bluewall geometry [13], analogous 
to the Schwarzschild tortoise coordinate r*. Parametrized thus, we see as in the dS 4 case that 
increasing B decreases the area, as long as we restrict the surface to lie within the future 
asymptotic universe /, i.e. /(r) > 0. As B 2 —* oo, these extremal surfaces become null with 
(u>() 2 = 1, corresponding to the past lightcone wedges of the boundary subregion, and have 
vanishing area. Thus extremal surfaces for a given subregion at X + can be constructed as in dS 4 
(Figured]) by joining two half-extremal surfaces: this is the blue wedge in region / in Figured] 
(the half-surface when not cut off continues as a null surface through the Cauchy horizon into 
region III , represented by the dotted extension of the blue line). As the subregion grows in 
size, this blue wedge approaches and eventually wraps the future Cauchy horizons. 

Figure 2: de Sitter “bluewall” Penrose diagram and 
some extremal surfaces with at least one end anchored 
at I + . The blue wedge is null, while the red timelike 
surface goes from I + to I~. 


asymptotic deSitter region 
T = 0 



asymptotic deSitter region 



One might imagine that there are timelike surfaces which are not restricted to just region / 
but instead start on Z + in / and cross over to II ending on the past boundary X~. These can 
be found with the parameter B 2 > 0 being finite. In this case, we see that ( w () 2 —* 0 as r —» 0 
and (ic() 2 y 1 as r y ~ near the horizon in I. Now after the surface crosses the future Cauchy 
horizon, we have /(r) <0 in IV. Requiring that (tc() 2 in (1341) satisfies (tc() 2 > 0 corresponding 
to real surfaces, it is possible to see ( e.g. by plotting as a function of r 0 r) that the parameter 
B 2 is bounded below by a critical value. There is a family of such surfaces: we will isolate one 
“critical” surface for a particular value of B , in what follows. Drawing analogies with the study 
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of the phase transition found in [03] (although the physical context there is different), we note 
that (w () 2 ~^ 00 when the denominator in (1341) approaches a double zero (with y = r 0 r), i.e. 



This corresponds to r c = -— ~ — which is just a little inside the Cauchy horizon in region IV. 


Note that {w'*) 2 \ Tc —> 00 here means this curve is normal to the w = const line here (these are 
straight spacelike lines passing through the bifurcation point and hitting the singularity in IV), 
or equivalently tangent to the r = const curve at r c in IV. The corresponding surface from 
r = 0 to r = t c can be drawn as a curve in the (r, w)-plane: it can be joined smoothly at r c 
with a corresponding curve from X~, resulting approximately in the red curve in Figure [2] This 
surface crosses the upper and lower Cauchy horizons at r = —, w = +00 and r = —,w = — 00 . 
The area of this surface has a leading divergence S ~ "of 77 “- Near the double zero, Aw acquires 
a large contribution and we can estimate S ~ -^V\Aw. This surface is vaguely reminiscent of 
the extremal surface in pHj which goes from one timelike boundary to the other: since the dS 
bluewall metric itself is related to the AdS-Schwarzschild black brane by flipping minus signs, 
it is perhaps not surprising that there exists a similar surface here (but timelike), albeit with 
no obvious corresponding interpretation. In light of ER=EPR [05], it is amusing to speculate 
that the subregion here corresponds to copies on both X ± possibly “entangled”, in some sense, 
thinking of the bluewall geometry as a “timelike wormhole” with the bifurcation region being 
the Einstein Rosen bridge. Note however that strictly speaking, all timelike geodesics go from 
Z~ to X + (unlike a shortcut in spacetime) either through the bifurcation region or through the 
Cauchy horizons, subject to the blueshift instability ra- 

With a w = const slice, real extremal surfaces likewise have 



(36) 


For B 2 —» 00 , these are again null extremal surfaces x 2 = 1 _^ r3 with vanishing area. These 
surfaces all lie on a w = const slice (thin black straight line from X + to X~ in Figure [2]). 

4 Discussion 

We have considered extremal surfaces in bulk de Sitter space (in the Poincare slicing) on con¬ 
stant boundary Euclidean time slices bounding subregions at future timelike infinity, motivated 
by the Ryu-Takayanagi prescription for entanglement entropy in AdS/CFT. Stemming from 
certain crucial sign differences, we have seen real extremal surfaces which are essentially the 
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restrictions of the boundaries of the past lightcone wedges of the subregion: these are null sur¬ 
faces with vanishing area. We have also seen complex extremal surfaces which do not always 
have real-valued area: this has parallels with analytically continuing from the Ryu-Takayanagi 
formula in AdS. In 0 LS 4 , the area is real-valued and negative. The area has structural resem¬ 
blance with entanglement entropy in a dual CFTd, with the leading divergence of the form 
Cd~n=r ■ the central charges Cd ~ i 1 ~ d 7 ^ s - here resemble those in the CFT energy-momentum 
tensor ( TT) correlators in [3] obtained in dS/CFT using Z C ft = 'k and a semiclassical approx¬ 
imation T ~ e lScl for the wavefunction of the universe. Alternatively, the strip entanglement 
entropy of a nonunitary CFT with these central charges would be of this form (assuming it 
exists), which the area of these complex extremal surfaces reproduces. This appears distinct 
from bulk entanglement entropy in de Sitter space (from a bulk density matrix via the wave- 
function T). We have also studied extremal surfaces in the dS black brane (where there is a 
finite cutoff-independent extensive piece), and the related dS bluewall spacetime. It is worth 
mentioning that there may exist other extrema of the area functional: for instance, we have 
required that x{r) parametrizing the strip width be real-valued, which suggests the path r = %T 
in complex r-space. This appears consistent with possible dS/CFT interpretations and also 
corroborates with our discussion of the dS black brane. However this may be restrictive and 
more general complex extremal surfaces may be relevant in complex r-space (see e.g. [53]). It 
may be interesting to understand if the analysis of [T8J can be applied in this case to obtain 
insights into extremal surfaces. 

While this analysis of bulk extremal surfaces could be regarded as simply a study of certain 
kinds of probes of asymptotically de Sitter spaces, it cannot pinpoint whether the corresponding 
area is expected to have a physical interpretation as entanglement entropy in the dual CFT, 
although the results do appear so, keeping in mind the central charges of the non-unitary CFTs 
in question. It is tempting to study this in light of the higher spin dS^/CFT^ duality of [4j. 
However the presence of massless higher spin fields might suggest that extremal surfaces which 
are geometric gravitational objects are not accurate (see e.g. mm which study entanglement 
entropy from Wilson lines in higher spin AdS holography). Nevertheless it is interesting to 
ask if these extremal surfaces have significance in some approximation where the higher spin 
symmetry is not exact. In this case, it would be interesting to explore the physical interpretation 
here more directly from a Euclidean CFT3 point of view. One way to think about entanglement 
entropy in held theory (lattice models) is in terms of the eigenvalues of a correlation matrix and 
a corresponding von Neumann entropy (see e.g. [T8] and more recently |49]). In that context, 
a simple model of a massless scalar held with wrong sign kinetic terms might suggest that the 
correlation matrix squared C 2 is related to that for an ordinary massless scalar held by a minus 
sign, so that C-eigenvalues A& become iX^. Then the associated von Neumann entropy is in 
general not real-valued: it would be interesting to understand this better. 
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